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In this paper we develop a method to describe perturbatively the entanglement entropy in a simple 
impurity model, the interacting resonant level model (IRLM), at low energy (i.e. in the strong cou- 
pling regime). We use integrability results for the Kondo model to describe the infra-red fixed point, 
ff) ■ conformal field theory techniques initially developped by Cardy and Calabrese and a quantization 

scheme that allows one to compute exactly Renyi entropies at arbitrary order in 1/T B in principle, 
even when the system size or the temperature is finite. We show that those universal quantities at 
arbitrary interaction parameter in the strong coupling regime are very well approximated by the 
_i ' same quantities in the free fermion system in the case of attractive Coulomb interaction, whereas a 

strong dependence on the interaction appears in the case of repulsive interaction. 
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I. INTRODUCTION 



Entanglement is a property allowed by quantum mechanics that describes the fact that generically, a quantum 
state of a system consisting of several subparts cannot be written as a product of states of the subparts. One way 
?h , to characterize, at zero temperature, this fascinating property of quantum systems, is to introduce the so-called 
entanglement entropy. Since it was first considered, this quantity has found many applications, ranging from the 
+J \ study of blakholes [1, 2], to modern methods of simulation of quantum systems [3], 

The study of quantum impurity problems from the point of view of entanglement entropy has led to many new 
insights and puzzles, both in [4] and out of equilibrium situations [5, 6]. In these problems, one typically considers 
a ID gapless bath coupled to a localized degree of freedom, and the questions of interest concern the entanglement 
of the two halves separated by the impurity [7-10] , the entanglement of the impurity with the bath [11, 12], or the 
entanglement of a region containing the impurity with the rest of the system [13]. In all those situations, one is 
Q ■ dealing with a bipartite entanglement entropy (associated to a partition of the system in two pieces A and B) that 
l— ~~ l [ characterizes to what extent the groundstate of the full system can be factorized into states for the subsystems A 
and B. Usually, the coupling of the impurity to the bath leads to a renormalization group flow, so the entanglement 
entropy exhibits different crossovers, and cannot be studied by the powerful methods of conformal field theory [14]. 
. While numerical results are available for several examples of the foregoing problems, analytical results are very few 
. Several difficulties are compounded here. One is the complexity of the entanglement entropy, which usually does 
not behave simply under the various Bethe ansatz and integrable quantum field theory tricks. Another, maybe less 
immediately obvious, comes from the fact the entanglement entropy, being a T — quantity, is very hard to access 
perturbatively because of Infra Red (IR) divergencies - its behavior in perturbation theory is in fact reminiscent of the 
■ screening cloud in the Kondo problem [15, 16]. Like for the Kondo model, the models we are interested in also involve 
flow to strong coupling at low energy: no matter how small the bare impurity/bath coupling 7 is, the asymptotically 
low energy states (below some scale T B (~f) that goes to zero when 7 — > 0) of the theory are qualitatively and deeply 
' affected by the coupling. In more technical, RG terms, the interaction between the impurity and the bath is a (maybe 
marginally) relevant operator, of scaling dimension D £ [0, 1]. The RG flow is characterized by an energy scale T B 
(that coincides with the Kondo temperature in the case of the Kondo model) , which marks the crossover between the 
weakly interacting regime, or ultraviolet (UV) regime, and the strong coupling regime, or infrared (IR) regime. By 
analogy with the usual thermodynamics quantities - such as the impurity entropy [ 7] - one expects the entanglement 
and Renyi entropies of the impurity, to bear a universal form and depend only on dimensionless parameters in the 
low energy (with respect to the cutoff - typically the bandwidth) limit : 

R imp _ Rn _ , ( v F v F 1 \ imp _ bulk _ / v F v F 1 \ 

N " R^- tN {LT B 'L T B ^ B '-) b ~ \LT B ' LqT b ' f3T B ' " J W 



A remarkable exception is the calculation in [I i] of the rate of entanglement entropy in the steady, out of equilibrium state of two wires 
connected by a resonant level, i.e. without interactions. 
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where L the length of subsystem A, Lq the size of the system, /3 _1 is the temperature, and where by "bulk" 
quantities we mean those of the system with no impurity 2 . Scaling analysis meanwhile shows that T B oc r y 1 ^ 1 ^ D . 
While the validity of the form (1) seems well established numerically [4, 12] at least in some cases, the natural 
attempt to calculate the universal functions perturbatively in 7 fails: indeed, while it is easy to set up the calculation 
and bring it down to the evaluation of explicit integrals, all orders turn out to be infinite [18], despite the presence of 
a natural IR cutoff - the length L of the subsystem. While these divergences simply mean that the scaling functions 
are not analytic in 7, they leave us facing the practical problem of evaluating a non perturbative quantity in a 
particularly difficult setup involving replicas etc. The present paper and its sequels are devoted to making progress 
on this important question, with the ultimate goal of computing entanglement rates in out of equilibrium, interacting 
situations. 

The impurity models we are going to consider have in common that the impurity degree of freedom is a two-state 
system (occupied or empty dot in the case of the IRLM, spin i for the Kondo model). At the IR fixed point, the 
impurity will end up being " absorbed by the bath" , the physical picture in a lattice regularization of the model being 
that it forms a singlet with the bath site neighboring the impurity [19]. In a continuum model, the hybridization 
of the impurity with the bath should be rather thought of as occurring in a region of size A B = v F /T B , the so-called 
"Kondo cloud" in the Kondo model. At zero temperature, one thus expects that the impurity has no effect on the 
entropies Rn,S when the partitioning of the system is done well beyond the impurity cloud (i.e. when the cloud is 
well inside region A, L ^S> As), resulting in the simple limits f N (0, 0, 0) = 1 and f s (0, 0, 0) = 0. Another simple limit 
at zero temperature is that of a disconnected impurity 7 — > + or T B — > 0, where one has 3 f N (oo, x, 0) = 2 1 ~ JV and 
/ 5 (oo, x, 0) = In 2. The entanglement entropy thus has the same limits in the UV and IR regime as does the boundary 
entropy s = \n(g) (in this general formula g is the "universal non- integer groundstate degeneracy") introduced by 
Affleck and Ludwig [20, 21] when varying the temperature between (3" 1 = 00 and = 0. In between those two 
simple limits lies a crossover regime where the interactions play at full and shape the crossover functions. It is probably 
tempting to speculate that, since they interpolate between the same values, the entanglement and boundary entropies 
could be the same functions, after formally replacing the length of region A by the inverse temperature. We shall see 
below that this is definitely not the case, and that entanglement and boundary entropies behave very differently, as 
a function of the interactions in particular. 

The Kondo situation is illustrated on figure 1 where we have represented the spin absorbed by the bath in the low 
energy (long distance) ("infrared" or IR) limit, and the spin decoupled from the bath at high energy (short distance) 
("ultraviolet" or UV). It is well known in these two, conformal limits, that the entanglement entropy of the region of 
length L with the rest of the system behaves as 4 

s = ^^ + s -l + lng (2 ) 

6 a 2 

(where a is a cutoff length scale, and s\ a non universal constant) so S lmp interpolates between In 2 and 0. 



2 Depending on the geometry, this system with no impurity may have a boundary. 

3 Note that the point 7 = is singular since then the impurity density matrix is totally arbitrary, not being specified by the contact with 
the bath. 

4 Recall also that the entanglement entropy of a region of length L sitting not at the edge but in the bulk of an infinite system takes the 

related form S = %ln — + si. 
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FIG. 1: Illustration of the low energy or long distance (IR) and high energy or short distance (UV) limits of the 
Kondo model, after reduction to the s-wave channel. In the low energy limit, the impurity spin is swallowed by the 
bath, while it is decoupled at high energy (denoted by a white circle). The black arrows denote right- and 
left-moving fields. A and B denote the partition of the system that is used to define the bipartite entanglement. 

The Kondo model (or the IRLM also considered below) being integrable, it is natural to wonder about determining 
the entanglement exactly using Bethe ansatz type techniques. This however remains a seemingly unreachable goal, 
despite progress such as in [22] . The next best option is to obtain what are usually extraordinarily good approximations 
using the form factors technique, directly in the continuum, scaling theory, after taking the appropriate "massless 
limit" of the usual massive construction [23]. This will in fact be discussed in our next paper. Meanwhile, there is one 
last tool we can try to put to use, following the calculation to lowest order performed in [12]: the expansion around the 
strong coupling fixed point. Such an expansion is bound to be better behaved as far as IR divergences are concerned, 
since one perturbs by irrelevant operators, whose correlation functions decay sufficiently fast at infinity to render all 
integrals convergent at large distance. Indeed, the lowest order calculation proposed in [12] relies on the effective 
Fermi Liquid theory for the vicinity of the Kondo strong coupling fixed point, and boils down, technically, to the 
calculation of a single integral of the one point function of the stress energy tensor (see below) , a perfectly well defined 
procedure, leading to a result that depends only on T B ,L as required. This procedure is however bound to fail in 
general, beyond the first order, since there are now strong UV divergences, and perturbation by an irrelevant operator 
does not determine a renormalizable field theory. However, in the integrable case, the existence of an infinity of 
conserved quantities does, in fact, provide a full control of the low energy hamiltonian: the necessarily infinite number 
of counter terms are all explicitly known [24], together with a well defined (analytical) regularization procedure. 
This means that, while we do not know how to calculate the scaling functions at small coupling (since they are not 
perturbative) , we have, in principle, a tool to determine them perturbatively at large coupling. Of course, only a few 
orders are technically manageable, but this is enough to gain an understanding of the scaling functions, and answer 
important qualitative questions. 

The paper is organized as follows: In section II, we introduce the impurity model and its representation in the 
ultraviolet and infrared limits. In section III, we present the method of infrared perturbation theory in the replicated 
theory needed to compute the entropies. We then turn to the actual calculation of the entropies in an infinite size 
system in section IV: our result for the entanglement entropy is summarized in Eqs. (42,43). Section V treats the 
case of a system of finite size 2Lq, our results being summarized in Eqs. (58,59). 

II. MODEL 

A. Interacting Resonant Level Model 

We start in fact not with the Kondo hamiltonian, but with the interacting resonant level model (IRLM). This is 
a simple impurity model that describes a tunnel junction between a localized resonant level at x = and two baths 
of free, spinless electrons. We note here that the model and all the manipulations below can be generalized to an 
arbitrary number of baths, but this will not change qualitatively the results very much. Later, we will discuss briefly 
the one- wire case. 

The IRLM also includes a Coulomb interaction (that is also called "excitonic" interaction) between the level and 
the baths. After a mode expansion in the baths and a linearization around the Fermi points, we end up with an 
hamiltonian H = Ho + H B , where H — J2 a =i 2 — * I^°oo dx^Pldxipa (we have unfolded the fermionic fields living on the 
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half-infinite lines into right-moving fields on the complete line) . The boundary interaction is given by 

H B = 7^(0)^+11.0. 

0=1,2 (3) 

+ U (V4Vi + 4^) (0)(dU - 1/2). 



ft is convenient to use spin-f/2 operators d) = j]dS + and $ d = S z + 1/2 to represent the impurity is a Majorana 
fermion). The tunneling anisotropy is parametrized with j^/2e ie ^ 2 = 71 + 172 and we consider even and odd fermions 
defined as a rotation of the original fermions 

( cosf sin§\ (i> x \ . , 

We take the bosonized form of this model with il) a {x) = i] a : e l ^' Va ^ . / v / 27r j with a = 1, 2, ±. Then (k q = %%) 

1 f 00 

#0 [^a] = 51 9 / ^ (^^a) 2 = #0 [</>±] (5) 

1 o * J— OO 



H B 

\/7T 



a=l,2 

7 



Efc (6) 
+ -±(d x <p+ + d xl p-)(p)s. 

where Ub = g{U) is a function of the capacitive term for fermions. This function depends on the regularizing scheme 
chosen to bosonize the initial model [25]. 

B. UV hamiltonian : anistropic Kondo model 

As in [26] one can absorb the interaction via a unitary transformation U —: e la ( ip++tf ' s, ^ Sz : with a = Ubl\f^ that 
relates the model to the anisotropic Kondo hamiltonian where the boson </>_ has decoupled 

H = H [<t>±] + : e lXrt >+^ : S+ + h.c. (7) 

where A 2 — — {Ub — it) 2 + 2tt and the bosons <f>± are again a rotation of the bosons ip± 

, ( 1-K-Uh -Ub \ / \ 

Hi #) (-) 

the scaling dimension of the tunneling operator is D — A 2 /8tt, thus for A 2 > 87r the tunneling is irrelevant and the 
model flows to its UV fixed point which consist of free right-moving bosons living on the complete line. 

In the following, we will be interested only in the interesting situation where the boundary term is relevant [D < 1), 
and deeply influences the low-energy physics of the system by driving it to a strong coupling fixed point: the model 
flows at low energy to the Kondo fixed point where the impurity hybridizes with the wires and we get again free 
bosons with different boundary conditions for + , namely c/>+(0 + ) = 0+(O _ ) + A/4, whereas we have trivial boundary 
conditions for 0_ , that is decoupled from the boundary. The renormalization group flow to strong coupling generates 
dynamically a new energy scale T B , that marks the crossover between the UV fixed point (7 = 0) and the IR fixed point. 
All physical quantities are universal functions of l//3T Bl l/LT Bl etc. This scale T B itself is of course non-universal, and 
it will depend on 7 and D. A simple scaling argument leads to the power-law dependence T B ~ 7 1 -- . 

C. The entanglement entropy 

The quantity we are initially interested in is the entanglement between a region (region A) of size 5 2L centered 
around the resonant level and the rest of the system (region B). Unfolding gives us two chiral wires connected by the 



5 The unfolding procedure doubles the spatial size of region A. 
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resonant level as illustrated on figure 2. After the unitary transformation, the entanglement can be decomposed into 
a sum of contributions from the cf> + and 0_ modes. Only the <j) + modes interact with the impurity as illustrated on 
figure 3. Finally, one may "refold" back the result to obtain exactly the Kondo geometry from figure 1. 



wire 1 



wire 2 



2L 



FIG. 2: The IRLM after unfolding the wires to get only right movers. 



o 
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FIG. 3: The <f>+ degrees of freedom correspond to a single chiral wire coupled to the impurity. 

Note that the geometry we consider does not correspond to studying the entanglement between the two wires 
connected by the resonant level [ ]. This situation is somewhat more complicated, in particular because the problem 
then cannot be turned into a purely chiral one - we will discuss it elsewhere [18]. 

D. Dual hamiltonian : approach of IR fixed point 

Since the boson 0_ decouples from the impurity, we are left with a single boson and in the following we simplify 
notations and set <f>+ = <p. Remarkably, as was shown by [27], due to the integrability of the anisotropic Kondo 
hamiltonian, one has a dual expression of the hamiltonian that is an expansion around the IR hamiltonian Hq R = 
\ dx (d x cj)) (altogether with its boundary conditions). The expansion takes the form of a series of irrelevant 
boundary operators 



H = H. 



n>0 ±B 



(9) 



The operators 2n + 2 are expressed only in term of a (holomorphic) deformed energy momentum tensor T = T — 
i\/2irad 2 (j), with T{z) — —2n : dip 2 : (z) and the coefficient a is related to the scaling dimension of the boundary 
perturbation via a = . The deformed theory with energy momentum tensor T(z) has a modified central charge 

The operators C?2n+2 are conserved quantities of the UV hamiltonian, they commute with 



l-6a 2 



1-6 



a-py 

D 



H and between one another. The first terms read: 

0- 2 = T 



e>4 = (t 2 ) 

O e = (0 2 4 ) 



12 



-(0 2 d 2 2 ) 



(10) 

(11) 

(12) 



We also know the value of the coefficients g 2 n+2, they read [27] 



_L_ (-1)" +1 (D/tt)" 
2^ (n + l/2)(n+ 1)! 



"r 


K 2(l-D)) 






r 


, 2(1-Z»))J 





2n+l 



r / 2n+l_\ 

p ( (2n+l)J 
V 2(1-1?) 



(13) 
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Note that considering the IR fixed point amounts to taking the limit T B —> oo, or equivalently to set all physical 
energy scales (temperature, inverse size of the box v F /L,...) to 0. 

We are now able to develop a theory of perturbations in 1/T B around the IR fixed point for the Renyi entropies. 



III. METHOD 



A. Renyi and entanglement entropies 

When a system is divided into two sub-systems A and B (more precisely when the Hilbert space can be written as 
J&A <8> J%b)> one can define the Renyi entropy of order N associated to the partition as: 

R N = R N B = Tr A ( PA f (14) 

where pa is the reduced density matrix of subsystem A obtained by tracing out B degrees of freedom, pa — Trsp 
and the partial traces Trx(-) — $Z|£)g ye x (£l ' l£) * a ' ie * ne t race on on ly a subpart X of the degrees of freedom. Note 
that Eq.(14) makes sense (i.e. the sum converges absolutely) for complex N with Re(iV) > 1 The Renyi entropy has 
the property R N B — Rfi A [28], hence the abreviation Rn used in (14). The bipartite entanglement entropy between 
A and B is defined as 

S = Sab = T?aPa In pa = Sba (15) 

One can deduce the entanglement entropy from the knowledge of the Renyi entropies by using S = — lim^r^! OnRn- 
The replica trick is a way to compute the logarithm : if we take integer powers of the reduced density matrix, the 
Renyi entropy can be interpreted (up to a factor) as the partition function of a big system consisting of N (coupled) 
copies of the original theory. Then, one analytically continues the result to arbitrary N £ C, Re(-ZV) > 1 to obtain 
iS (see Ref. 30 for considerations about the non-existence of ambiguities in the continuation - essentially because 
TV = oo is an accumulation point). Note that the replica trick 



B. The replica trick 

We now specify the partition A/B, and choose a spatial division of the system into a single finite interval and its 
complement, i.e. xeAii—L<x<L and x <G B otherwise. In the case where the initial full system has conformal 
symmetry, it was shown in [29, 30] that the replicated theory bears a very elegant structure and can be described by 
a In orbifold theory [31, 32]. 

To arrive at this, Calabrese and Cardy first work with an Euclidian path integral representation of the full density 
matrix p at temperature If the action depends on a collection of fields : S [{<p(x,r)}], a matrix element of the 

density matrix p is written: 

which is represented graphically by figure 4a. Note that in our case, contrarily to [29, 30], the action S = So + S B is 
not the free CFT action Sq. 

A matrix element of the reduced density matrix pa is labelled by two spatial configurations of the field <^ 2 ) ( x ^ ^) 
which live only in the region A : ((ffi IpaI^i /• Since the partial trace amounts to take (f>(x eB,0)= <j){x € B 7 (3) and 
to sum over the fields <p(x £ B), graphically this is represented by a gluing of the boundaries at r = and r = f3 in 
region B, leaving a surface with a cut along x € [-L, L], over the lips of which the fields (f>f^(x) live. This is shown 
at finite and zero temperature in figure 4b and figure 4c. 

Performing the last trace over A yields a compact expression for the Renyi entropies Rn that can be expressed as 
a ratio of partition functions: 

Rn = (If (17) 

Z N = I V<j>e~ s W. (18) 

J-Rn(P) 
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— ^ t = /3 

02 

(a) A matrix element of the full density 
matrix p 



B r = 



B 



- x 



(b) Reduced density matrix at finite temperature /3 1 . 



B 



,r T = OO 

(c) At zero temperature /3 = oo. 



FIG. 4: Graphical representation of full (figure 4a) and reduced (figures 4b and 4c) density matrix elements 
((f>i(x)\ p(A) \4>2{x)). The fields <fn, 4>i that label the matrix element are the boundaries in the path integral formalism. In this 
formulation, the partial trace over B amounts to gluing together the boundaries in region x £ B. The fields <f> propagate with 

the complete action S. 




FIG. 5: The TZn geometry. The thick line at x — is the integration contours of the boundary operators we insert. In this 
example we have inserted two operators, coming from two different copies of the system. 



Here Zn is the partition function of the theory living on the surface TZn(P) that is obtained by gluing cyclically the 
N open cuts of figures 4b and 4c, so that 1Zn(/3) is the ./V-sheeted Riemann surface TZn when (3 = oo, and a more 
complicated surface of higher genius at finite temperature. 

The TZn geometry We now see the geometry in which we will have to compute correlators of boundary operators 
in order to get the Renyi entropies. The matrix p% is a fusion of the cylinders/planes via the "lips" of the cuts, 
and taking the trace over region B amounts to fuse the last one with the first one. This gives the TZn geometry, 
shown in figure 5. In the path integral formalism this means that we take N copies of the system, we have a 
collection of fields (f)(j\ , j = 1...N, and two inverted branching points, one at z = —iL and one at z = iL, so that 
(/>(j)(0 + — ix) — (f)/j_i)(0~ — ix) for —L < x < L. 



C. Infra-red perturbation 



1. Inserting operators on TZn 



We want to make a perturbative expansion around the free action So : S = So + S B with S B = J dT_ff B (0,r) and 
we develop in the path integral (16) the exponential 

exp(-S B ) = £ J dn...J dT n ^J-Y[H B (0,Ti). (19) 
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Graphically a matrix element of p is represented within this expansion in figure 6. 



01 




X 


Hb(oA) 


T 






<t>2 



T = 



FIG. 6: A matrix element of the full density matrix p, here at order n = 3. The fields cj>t, <j>2 that label the matrix element 
are the boundaries in the path integral formalism. Propagation of the fields <f) with the free action So and insertion of 

boundary operators. 



r imp 

'A — — ~ R NT^—Tn ( 21 ) 



Expanding these path integrals around Ho, one gets an expansion in powers of 1/T B of the Renyi entropy and thus 
of the entanglement entropy. The expansion is of the form 

R N = j XV + XV-e~ s / ( [ Vip+Vip^e-'A (20) 

J-Rn \J7li J 

<e- s «) Ki * N (Z™v)" 

where i?™ is the unperturbed Renyi entropy, and is obtained from Eq. (20) by replacing the action S by So 4>-] 
(i.e. setting H B — 0). 

The expectation value (-) Kjv is taken in the theory living on the Riemann surface TZn- 

(X) Un = j ^ + le" SoI * +1 / f V<j} + e- So[4, + ] (22) 

N JTLn JTIn 

where we have factorized the odd modes since the interaction H B does not depend on 0_ . 

Our computation is then a standard perturbative expansion, but it is carried on a complicated worldsheet. There 
are three factors in Eq. (21): 

• The quantity is the Renyi entropy at the conformal point. It has been computed in [30] and we will come 
back to it in the next section and in the appendix A. 

R% = c N (^-J (23) 

5 IR = - lim d N (J$) = °- In ( 2 -f) + d x (24) 

with c = 1 the central charge and cn and are non universal constants. 

• The denominator in (21) is nothing but the impurity partition function of the system raised to the power 
N, where the impurity partition function is defined as the ratio of the partition function of the full model 
Z = r £re~ l3{ + H+H ^ by that of the system as the IR fixed point: 

Z 



2^imp 



Its computation involves a series of expectation values on a single-sheet geometry. For example, in the case of 
infinite size and zero temperature, the expectation values are taken on the complex plane and we will see shortly 
that the denominator is exactly equal to one. 

• Similarly, the numerator can be written as a ration of partition functions, but for the system living on a A^-sheet 
geometry: 

Z T^^-(^) nN (26) 
We will see in the next section how to compute it perturbatively. 
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2. Explicit expansion 



We now turn to the description of the pcrturbativc expansion of the Renyi entropies, or equivalently the partition 
function on the Riemann surface 1Z jy of the full theory with boundary interaction (9). 

The partition function Zjy is now an infinite sum of contributions involving fc-uple integrals (see figure 6) of fc— points 
correlators on TZn of operators O ni , (j\i = i O ni (zi)^, where Zi — Tj is a point on TZn at position x — 0. To start with, 

the exponential e~ Ss is time-ordered in imaginary time. 

Recall the effective IR hamiltonian is obtained from the expression of the boundary state describing the full interact- 
ing theory in the closed-chanel geometry [27], whose expression is naturally regularized by inserting the perturbating 
operators at points with < —x\ < —X2 < ... < — Xk <C 1. Note that this regularization of the UV 

divergences by point splitting is the only one consistent with integrability. Within the point splitting regularization, 
the fc-uple integral is done over paths that are slightly displaced and therefore never cross. 

Furthermore, recalling that the perturbing operators O n are conserved quantities, one observes that we can freely 
exchange the order of two contours. Indeed the term that is generated when two contours i and j are exchanged is 
simply ... J dzj (fi z dziO ni (zi)O nj (zj)... = (here ... indicates that this expression appears inside some correlators and 
some integrals). The vanishing of this expression resulting from the 0„'s being commuting conserved quantities, that 
translates in the OPE: O n (z) ■ O m (uj) ~ ... + dO n>m {ui) + reg. with 0„ jm a local operator. 

We are now in a position to write down the expansion of the impurity Renyi entropy -§nf : 

^ P -S = E C fM^ {M (27) 

1 ' rt 



l=1 J{x=0} \ i 

C {n t } _ ^ jjfc j.-E?=l( 2 "i+ 1 ) ~Q (g2m,+2)' ?i ^ 




qii 

where the sum in (27) is taken over all integers vectors n = (0 < n\ < .... < nk), and one has introduced the integers 
mi,qi writing n = {mi, mi, m-z, m2, ••• m p , m p ) with m, < m i+1 . For fixed n i; the integral over z^ is to be 

91 92 q p 

taken along the N possible disconnected paths on IZn going from — oo to +00 in imaginary time, and position x = 
(thick line in figure 5). 

In the next section we will see how we can compute parts of the A^ ni ^ in the case of infinite size and T = 0, where 
the sheets of IZn are the complex plane. 



IV. INFINITE SIZE 



In this section we report results for the perturbative expansion in the simplest case of a system of infinite size and 
at zero temperature. We use an orbifold representation of the expectation values that appear in (21) and conformal 

transformations to explicitly compute the expansion up to order ^ ^ 6 . We refer to appendices A and B for details 

about twist fields and orbifold representation. 



A. Infrared fixed point 

We start by briefly describing the theory at its IR fixed point (formally obtained by setting T B = 00) focusing on its 
equivalence with a Zjy orbifold theory. This method was introduced in Ref. [30] by Calabrese and Cardy to compute 
the first term . It relies on the introduction of "twist fields" which are inserted in a path integral on a more 
simple geometry and which encode the twisted geometry oUZn- 

Comparing (Al) and (22), Calabrese and Cardy find that j nN Vip + e~ s ° oc ($(u)$(u)), where are the twist 
fields, that are inserted at points u — —iL and v = +iL. This leads to the well-known result for the Renyi entropy of 
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an interval of length 2L with the rest of the system in a free theory 

2 £\ -vi N ~w) 



R% = C "[—J (30) 

where we recall that we have chiral fields but two copies <j>±, resulting in the factor c/6 in the exponent. When taking 
the limit N — > 1, this results in the by- now well-known expression [2, 29, 30] for the entanglement entropy at the 
strong coupling fixed point: 

S m = — In ( — ) + non universal (31) 



3 V a o J 

with c = 1 in our case. 



B. Away from the IR fixed point 



A generic term in the perturbative expansion corresponding to the insertion of k boundary operators, consists in 
a fc-fold integral on the surface TZn- One can actually deform the paths of integration continuously, one by one: 
the integrands are k— points correlation functions of the Oi's, therefore analytical everywhere on (TZN) k , except for 
the branching points at z = u, v = —iL, +iL. Bringing the contours past the singular point z = v leaves behind a 
contribution that is simply a path Pn circling ant i- clockwise the singular point z = v on TZn ■ The other contribution, 
being past z = v, can be moved to infinity, where it vanishes due to the clustering properties of the fc-point correlator. 
This is shown in figure 7 in the case N = 3 and one contour (if there are many contours, one obtains nested integrals 
around z = v). 




FIG. 7: Deforming one contour on TZn, illustrated with N = 3. 



In fact, since the path of integration Pn is invariant under cyclic permutations of the replicas, we are dealing here 
only with very special operators, namely Z^r invariant operators X^(z) = Ylj—i X(z^) = ^2jLiXj(z), with j a 
sheet index. The equivalence between the theory defined on TZn and an orbifold theory allows to express the mean 
value of any such operator as: 

/ rnw (X mh (z)<i>(u)${v)) , x 

where expectation values in the right hand side are taken in the orbifold theory on the plane. Hence, introducing the 
notation Res [A orb • $] = {X olh ■ $} x for the terms in ^ in the OPE A orb (z) • $(w) = £ n£Z * w {A orb ■ 
one arrives at the point where one can represent, at the operator level, the integral over the path Pj\r: 



$(w)Res[0° rb • $] (v) 
dz(O n (z)) nN =2m^ '- (33) 



This is easily iterated (the different contours are nested) with the result: 



r , , Mu) Res n $ (vY 
A^(u-v) = (2in) k + — ^ — '- (34) 



Res s [$] = Res [0 2ni +2 ■ - Res [0 2nk+2 •#]...] 

(35) 
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Considering for a brief moment terms with only O2 insertions (i.e. n\ — TI2 — ... = fik = 0), and 
using the fact that Res [T or b • -X"] = dX for any local operator X, one immediately finds A^° °^(u — 

v) = /*(u)*(t>)^~ (2i7r) fc 9,^$(u)$(w)^ so that Efc>o ( ~ g2 /, TB) " A^ - } = /*(u)$(t;)^~ - Sfsaa )\ = 

I 1 + J where one has used that gi — — \, u = ~iL and v = +iL. One can thus resum the full series 

for those terms consisting only on C 2 contributions. 

More generally, since an insertion of O2 amounts to a derivative with respect to the point v, one can resum all O2 

subdiagrams, J2 k >o ^ {ni '-'" r '°"' 0} (u-«) = + 2i/T B )\ x A*™ 1 -- (it - u - 2i/T B ): 

This allows for a partial resummation of the perturbative expansion, all O2 insertions simply amounting to a shift 
v — > v' = v + 2i/T B . We thus arrive at the final compact form: 

R "" = ( 1 + ik)"" *»(•— I 



5 imp = s 0a + s ( u _ „ _ |l ) (36) 



= iln(l + z L) (37) 

where one has introduced the reduced impurity Renyi entropy free of O2 insertions: 

R N {u-v') = C {ni} A^(u-v') (38) 

n \ iii>Q 

and S = — limjv-yi OnRn- This rewriting considerably reduces the number of diagrams that need to be computed to 
achieve a given order in the expansion: only diagrams free of O2 insertions have to be considered, and at the end of 
the calculation, one replaces v' by v' = v + 2i/T B = iL + 2i/T B . 



C. Results 



We now turn to the actual evaluation of the reduced partition function Zn(u ~ v). According to Eq. (35), one way 
would be to compute nested residues of the perturbing operators O4, Oq, .... The situation is not as simple as for 02- 
For instance if we want to compute in the orbifold theory an expectation value such as ^ ((T 2 )(zW)) TC = 

Si ((-^; 2 )( z ))-r, the orbifold version of ^ (T?) is not (T 2 rb ) and in particular we have a priori no clue of the OPE 
of the former with twist fields. 

It is nevertheless possible to compute such an OPE (see appendix C for the explicit form of Res [C° rb • <£>] ) , but 
it notably involves to know in advance the correlators we are at the end interested in. Contrarily to the effect of 
O2 insertions, there is (probably) no simple differential operator representing higher operators O n at the level of 
correlation functions, presumably forbidding a possible resummation of the Renyi entropy in the same spirit as could 
be done for C 2 contributions. 

The contributions A^ ni>0 ^ can nevertheless be evaluated in a straightforward way: any such contribution can be 
obtained from the knowledge of the n— point function of the operator O2 on TZn, by then performing the relevant 
normal order operations on the surface IZn (this is the route we have chosen, see for instance appendix Bib). Another 
equivalent way would be to derive directly the form of the orbifold operators on TZn- For instance, the expression of 
the operator Of b (z)= £\ 4 {z { ' } ) can be deduced from the transformation of the A = (T 2 ) — fjj<9 2 T quasi-primary 
operator, that reads 



Kn N ( Z ) = /» 4 A c (o,) + . C ALA (y (z) 2 TcM + ^f) (39) 

upon the conformal mapping lj = f(z), z € 1Zn,u £ C that uniformizes the surface TZn into the complex plane. 

We now turn to the results up to order 7 in l/(LT B ). The partial resummation (37) of O2 insertions leads to a signif- 
icant simplification of diagrams that need to be computed, from 17 diagrams 

/if 1 . 1 } 4(0,0,0} /|{0,0,1} /i{0,0,2} ^{0,1,1} a{0, 0,0,0} a{0, 0,0,1} ^{0, 0,0, 0,0} ^{0,0,0,0,1} a{0, 0,0, 0,0,0} a{0, 0,0, 0,0, 0,0} 
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to only 3 reduced diagrams (i.e. diagrams free of O2 insertions) and A^ 1,1 ^ at order 7. After tedious 

calculations, we find: 

Rn(u-v) = -^A^-(u-v)-S'A^ { u-v) + M r A^Hu-v) + 0((LT,)-') 



A^{u-v) 



-inlN 2 - 1 



24(2 - D)(l - 2D) ((1 - L>) 2 (334iV 2 - 246) - D{64N 2 - 31)) 



5760D 2 7V 6 (u-i 
+LW (-18L>(1123iV 2 - 667) + (1 - D) 2 (7iV 4 1193(W 2 - 5697)) 



A^(i 



(N 2 - 1) n 2 
N 6 (u - v) 6 



(1~ADY{A~D) 2 / . „ .3 , 144 



5(7^ _ \f + jv(-36 + 79(AT - 1)) 



2880L> 2 V '13 



, 3521D 2 - 1576D (l + D 2 ) + 216 (l + D 4 ) 
+9iV(20- 29iV 2 + 13Af 4 ) \ J \ ) 



(40) 



1820 D 2 

Due to the prefactors (iV 2 — l) 2 in all but one of the contributions to Rn, the reduced entanglement entropy bears 
only one contribution up to order 7, that can be expressed as a function of parameter a = 

S(W -*/) = ! ( 4 « 4 - ^ + 9) + O ((LT B )-») (41) 

where u',v' are the generic points where the twist fields are inserted in order to compute the quantities A^ 71 ^ with 
n.i > 0. Now, taking u' = —iL, v' = iL + 2i/T B in this expression and keeping only terms that are O ((LT B ) -7 ) leads 
to the expression for the full entanglement entropy: 

S = S m + S imp (42) 

s '-> " i >" i 1 + ik) " i <5H5F (' " S;) (4q4 - 8a2 + 9) + 0(<i^ ■ ^8, (43) 

with 5 IR the entanglement entropy at the IR fixed point given by (31), and we recall the I?— dependence of the 
coefficients 54 and a: 

_ D ( T(D/2(l-D)) \ 3 r(3/2(l-£>)) (1 - D) 

94 6tt 2 V T(l/2(1 - D)) J r(3D/2(l -£>))' ' 1 ' 

Note that in (43), the first term in the right hand side has to be truncated at order 7. 

The second term in the right-hand side of equation (43) is the first term in the perturbative expansion that depends 
on the strength of interactions U in the original model. We note that its sign is negative, and this is consistent with 
the fact that the impurity entanglement entropy should eventually saturate at the value ln(2) in the limit LT B — > 0. 
An interesting observation is that this correction is minimal at the free point U = in the IRLM or the Toulouse 
point in the Kondo model, see figure 8. 

In the regime of attractive coulombic interaction, i.e. at negative U corresponding to D > 1/2, the correction (41) 
displays a very weak dependence on D, and the entanglement entropy, at order (LT B )~ 7 , is essentially that of free 
fermions {7 = 0. 

On the other hand, at very small D (this situation can be reached in the one-wire IRLM by choosing a sufficiently 
large repulsive coulombic interaction U), the correction (41) scales as 1/Z? 4 and is thus unbounded: the universal 
crossover function defining the entanglement entropy is thus deeply affected by the interactions for U > 0, even in 
the strong coupling regime. 

It is instructive to compare the universal crossover function for the entanglement entropy with that of the boundary 
entropy s introduced by Affleck and Ludwig [20, 21]. In the UV limit, both those quantities reach the value ln(2). 
In the IR limit, they both vanish. The boundary entropy s is simply related to the partition functions of the 
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FIG. 8: Impurity entanglement entropy: first correction to the contribution | ln(l + 1/LT B ) stemming from 02, displayed 
here at fixed LT B — 1 as a function of the interactions. D = ~ corresponds to the free point, where this correction is minimal 
(see inset). A significant enhancement of this negative correction is visible when reducing D. 



system at finite temperature (3 1 : s = (1 — 0dp)hiZ imp with Z lmp = ^of^ an d the partition function of 

the system without impurity. One immediately gets lnZ lmp = \a{eT Ss ) „ where the expectation value is taken 

at temperature i.e. on a cylinder of circumference ft. Because the perturbing operators O n are commuting 

conserved quantities, in the perturbative expansion of this quantity around the strong coupling fixed point, the 
contours can be deformed and moved one by one to infinity along the cylinder, resulting in (e _SB )^ = e~' Sl1 '' 3 and 

finally lnZ im P = - (S H } = -0 £„ ff2n+2 T B - (2n+1) (O 2n+2)/3- The boundary entropy thus has an expansion in odd 
powers of T/T B . Evaluating the one-point functions (027i+2) ( g in the cylinder geometry with circumference /3, one 
arrives at the expansion 

1 / 7T \ 9-4a 2 ( 7T \ 3 96a 4 - 340a 2 + 425 / tt \ 5 _ //1/omm 

s = -^ 2 3 [ws ) ~ m ) ~ nm — ^ — (wJ + ° ((1//3Tb) } • (45) 

We note that a dependence in the interaction parameter D appears at order 3 already, making the boundary entropy 
more sensitive to interactions than the entanglement entropy. 



V. FINITE SIZE 



We now want to compute perturbatively the expression of the entanglement entropy around the IR fixed point but 
in the case when the size of the system is finite, of length Lq in the original geometry or 2Lq in the unfolded geometry. 
In this case, the replica trick leads to expressions of Renyi entropy similar to Eq. (20), but instead of considering a 
surface TZjq made with N infinite planes with branching cuts, one gets a surface Cn made of N cylinders periodic 
along the x direction with the same branching cuts between them. 




FIG. 9: The reduced density matrix is now a cylinder, periodic in space, with a cut along the x axis. The Cn geometry now 
consists in N cylinders connected along the cuts. The partial trace over subsystem B results in the gluing of two semi-infinite 
cylinders along the dashed line, at r = 0. The impurity is inserted at x — (not represented). 



14 








\ H B {osi) H b (uj 2 ) 




'* 

X 


x X 

A(«) 



FIG. 10: The orbifold theory is now living on a cylinder and one has to uniformize it to the complex plane. The integration 
contour corresponding to the impurity insertion is then radial and goes from uj — to uj — oo, passing between the twist and 
anti-twist fields. The dashed line corresponds to the partial trace over subsystem B. 



A. Infrared fixed point and twist fields 



The first contribution to the entanglement entropy at finite size in the strong coupling limit LT B — > oo is, as before, 
the entanglement entropy at the IR fixed point. This quantity, that was computed by Calabrese and Cardy [30], can 
be obtained as the two-point function of the twist and anti-twist operator at finite size ; since those operators are 
primary operators, by using the conformal mapping that sends the cylinder onto the plane: 

z e C — > uj = e^ z e C (46) 

one easily computes ($(tt)$(«) ) , yielding: 

S^lnf^sin^W (47) 
3 \na Lo J 

Here ao is the inverse of the ultraviolet cut-off (e.g. in a lattice regularization, it coincides with the lattice constant), 
and c[ is again a non universal constant. 

The finite size infrared fixed point is now our reference point around which one wishes to expand the Renyi and 
entanglement entropies. The orbifold representation allows one to compute the expansion even when the surface is 
Cn- Indeed, since this representation gives the correlators of the original A-sheeted geometry in terms of correlators of 
the orbifold theory in the geometry of a lone sheet, one only has to perform the change of variables from the cylinder 
to the complex plane (46). 

However, there is a slight difference with before. When considering the surface TZn, one can move the integration 
contours freely in the orbifold representation and only get a contribution from residues around one of the u, v points, 
and this fact greatly simplifies the computations. In the finite size case, we first have to flatten the cylinder and 
change the integration contour accordingly (see Fig. 10). In the new geometry the integration contour on the plane 
is not a circle anymore, it is instead not not closed leading to more complicated integrals. It results that the insertion 
of the operator Oi has no simple action and we cannot derive an equality such as Eq. (37). 



B. IR divergences 

At finite size, one also has to be careful to take into account the expansion stemming from the impurity partition 
function Z lmp — ( e ~ S ") Kl (the denominator in Eq. (20)), which is now not equal to 1 because of the non trivial 
expectation value of O2, and of higher order perturbing operators, on a cylinder. 

This additional complication is easily understood by just considering the first order perturbing operator O2: having 
put the system on a cylinder, the curvature induces a non- vanishing local energy density, (O2) = (T) = eo = — #75. 
It results that the impurity partition function is now formally infinite : introducing a (small) finite temperature /3 _1 

for a moment to regularize those divergences, one has Z lmp = lim^oo e~ eoTr o , so that each term of the expansion of 
Z lmp in inverse powers of T B is formally infinite at vanishing temperature. 

These divergences will also appear in the A— sheeted geometry TZn, when computing -Z™ p . The impurity partition 
functions Z lmp and Z™ p cannot be defined on their own and the quantities A^ ni ^ defined in (28) as constituents of 
Rn are now ill-defined. 
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In fact these divergences compensate term by term in the expansion in 1/LT B of Rn / R™'- the ratio between the 
impurity partition functions is perfectly well defined yielding a finite answer for the Renyi and entanglement entropies. 
One way to see this is to take the logarithm of the Renyi entropy : 

lnR N - lni2$ =ln(e- s *) Cjv - Nln(e' s -) Ci = {e"**)^ - 1 - ^((e"^ - 1) (48) 

where we used the well-known relationship between the generating functionals of connected and disconnected expec- 
tation values of operators. The connected expectation value of a product of operators O ni (zi) is defined recursively 
as 

(O n (z)) c CN = (O n (z)) CN (49) 

and 

(o ni (zi)... o nk { Zk )) CN = (o ni («!)... o nk ( Zk )) c CN + j2 < n °^ ( z i)) c CN < n °»> fe)>L ■ • ■ ( 5 °) 

part. i£a j£/3 

where a, /3, . . . are non trivial subdivisions of a partition of [l,k]. 

This writing allows to regroup the terms of the expansion in 1/T B in well defined quantities. Indeed, 
these connected (or non factorizable from the point of view of the integration variables) contributions of 
iO^{zi) . . . $(7i)$(u)) c /($(w)$(w)) c can be written as the sum of two terms. 

The first term depends on u, v has to vanish when u — > v (or u, v — > ^iLo, cutting open the cylinders) and encodes 
the local influence of the insertion of twist fields at u,v. This first contribution does not produce infinities when 
sending all the Zj — > oo (while keeping u and v finite) because of the clustering properties of the connected correlator 
between the O ni and $, $. 

The second contribution which does not depend on the points u, v (from (O^izx) ■ ■ ■ } c ) can be shown to be 
divergent when all the points are sent to infinity. In this limit, the fields inserted at points Zi do not feel any 
longer the cuts on the cylinders, and the correlator is asymptotically that of fields living on N disconnected cylinders, 
yielding a factor N, the number of copies, for the result of the integral. 

In a nutshell, it can be shown that \(~SB) k ) r = (u — v)Fi(u,v) + N x diverging part, where the diverging part 

is independent of N and is exactly compensated by N((—SB) k ) c ■ This shows that the expansion in lni?jv — lni2™ 
are well defined. 

As a result, when computing the Renyi entropy Rn at finite size, one only has to take into account the connected 
correlators in Z™ lp , and do not care about the diverging parts that are in fact all canceled out by the denominator 
2 lrap at every order in the development in 1/T B . 

C. Twist fields at finite size 

We can now easily generalize the calculation of the first order terms to the finite size case. Under the conformal 
mapping (46), the singular points u,v of the are mapped onto v — e J7ri / Lo and u = e~ l7rL / L °, so that one has for 
instance: 

^ X >c N ($(u)$(«)) c 



i = l 



((£) T olh (u>Mu)*(v)) c c (dz\ 



-2 



- 



{z-uj} 



The last term yields a diverging contribution that is exactly canceled out by a term coming from (z imp ) N , the 
denominator of Eq. (20), so that we get: 



7T 



L o Jo 



hbj{u-v) 2 N 



{lo - uf (u - 5) 24w 



^h{L/L ) + div. (52) 



with the function /i(a) = a(l + 7r(l — a) cot7ra). Similarly one can obtain the higher order contributions involving 
insertions of the perturbing operator O2 (see Appendix B 2). 
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The next step consists in computing correlators of higher order perturbing operators C„>2, on the surface 
Cpf. One route would be to express the Z^-neutral combinations corresponding to those operators, namely 
G° rb = Ylj=x ®n ( z ^)> m terms of orbifold (quasi-) primary operators. We have chosen another equivalent route, 
namely a quantization method that allows to directly compute correlators of more complicated operators like O4. 

D. Direct quantization 

A direct way to compute the correlators Ay 1 *' at finite size is to look directly at the structure of the Hilbert space 
in the twisted sector and to consider a quantized version of the fields. This is a generalization of the quantization of 
fields with Neveu-Schwarz or Ramond conditions. The field living on N sheets is in fact a N— valuated field and we 
shall quantized it this way. 

1. Monodromy relations 

To be general, in the replica trick we consider a CFT on a Riemann surface with N sheets. Each sheet is described 
by the complex coordinate z = r — ix. We consider a real free scalar field ip(z, z) living on the Riemann surface, its 
value on a given sheet j = 1, . . . , N is denoted by ip^ (z), so that we have in fact N copies of a (chiral) CFT. 

These copies are not independent : we consider a branch point at z = 0, such that (p(j\(0 + — ix) = <yS(j_i)(0~ — ix), 
for all x > 0. In other words, when one turns anticlockwise around the point z = 0, one decreases the index of the 
sheet by one : 

m {ze 2i *) = y {j - l) {z) (53) 

making afterwards the conformal transformation (D5) will allow one to consider the branch cut that we are interested 
in. It is more convenient to consider the Fourier-conjugate fields, in order to diagonalize the monodromy relations. 
We thus introduce the fields tp^ (z, z) = J2f=i e 2 ** 7rjfe / JV <P(j) (z, z), which verify the relations 

<p( k )(ze 2in ) =e 2 ^V fc )(z), (54) 
and since the holomorphic energy-momentum tensor is given by T(z) = — 271"^^ : dtp>a\dip(j\ : (z), we also have 

T(z) = -2tt^ : d V ( - k) d V (k) : (z) =^tW(z). (55) 

k k 

Note that <^~ fe) = ip^ N ~ k ^ = ip^* : the index is to be understood modulo N. 6 



T 




FIG. 11: Fields living on a N— sheeted Riemann surface. The cut is along the x axis, between u — 0, v = 00. 



If N is even, we now have a free theory with two real bosonic fields (<p(p) , Wjv/2) an< f the associated energy-momentum tensor is 
rp(o) ^rp(N/2)~j an( j (j\T — 2)/2 complex bosonic fields (the associated energy-momentum tensors are the +T( -fc ' = —4Tv(df^df^~ k ^)). 
One the other hand, if N is odd, we have only one real field <P(o) an< i (N — f)/2 complex fields. 
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2. Mode decomposition 



A canonical quantization on a cylinder altogether with imposing the boundary conditions (54), leads to the mode 
decomposition 



^ -o®*-™ m € (Z - k/N) 

V47T ^ m 



with the commutation relations 



k+k' Om-\-n 



(56) 



(57) 



Acting on the twisted vacuum we have a™ | a) = if m > 0. Here |cr) = (^^Lq 1 |<J^ fc ^) is the twisted vacuum and 
the Jcr^) are the highest- weight vectors of the representation of the algebra (57) ; k is the index of the monodromy 
sector. 

The twist field $(2) is the primary field associated with \a), and it can be written as a product of the twist fields 
associated with |<jW) : $(z) = n' fc $W (2). The prime on the product simply denotes that we count k only for the 

independent fields. We shall see in the next paragraph that the conformal dimension of is hk — (2) fc ^^ , the 
factor 2 depending on whether k is the index of a real or complex field. In the twisted sector the vacuum expectation 
values are abbreviated as (X) = ^tht^-. 

. Y" a > 

Within this canonical quantization scheme, one can build the perturbing operators O n and evaluated arbitrary 
correlators thereof by using Wick's theorem for the bosons <p( k ' . Finite size is eventually taken into account by means 
of the conformal mapping (46). We report on the details of the actual computations in Appendix D and directly turn 
to the end result for the entanglement entropy. 



E. Results 

Finally let us summarize the results for the entanglement entropy of a segment of width 2L in a system at finite 
size 2Lq, up to order 0(T~ 4 ). The net expression for the entanglement entropy is: 

S = Sf +S imp (58) 
' h{L/L )- — 1— / 2 (L/L ) 



6LT B 1 ' 12(£T t 



V.) 



1 ' ^ ^V3(Vio) + ^7^(9-4a 2 )(L/L ) 2 /i(iAo) j (■><>> 



\2LT B J V9 120 

where S™ is the Calabrese-Cardy entanglement entropy (47) in the IR limit LT B — > 00, and the impurity entanglement 
entropy is expressed in terms of the functions (we also give their asymptotics when a = L/Lq —> 0): 

fi(a) = a(l +7r(l - a) cot att ) - 1 - 7r 2 a 2 /3 + ttV/3 + C(a 4 ) (60) 
/2(a) = /1(a) 2 + a 2 (l - a)V - 1 + 7r 2 a 2 /3 - 4ir 2 a 3 /3 + C(a 4 ) (61) 
/ 3 (a) = /i(a)/ 2 (o) + 2a 3 (l - a) 2 7T 2 ~ 1 + 7r 2 a 3 + G(a 4 ) (62) 

We can see that the interactions (which are hidden in the coefficients a — and 54) do not affect the entanglement 

entropy before the third order in 1/LT B : the effect of the interactions at strong coupling is thus a subtle one, just 
as in the infinite size case. We again note that this conclusion is invalidated at strong repulsive Coulomb interaction 
in the one-wire IRLM, where D can be made as small as wished: in this limit, the term proportional to 34 scales as 
1/D 2 and thus eventually becomes dominant at sufficiently strong U > 0. 

Our result for the entanglement entropy has the expected property that l S lmp vanishes order by order in the limit 
L — > Lq: in this limit, region B shrinks to a point (and this limit is not singular since region B does not contain the 
impurity) so that all Renyi entropies evaluate to 1. 



VI. CONCLUSIONS 



The entanglement entropy is an elusive quantity, which in many respects differs from usual local observables. In 
1+1 dimensional many-body systems at their critical point described by a CFT, it has been shown that the replica 
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method allows for a computation of the Renyi and the entanglement entropies. In this paper, we have shown how 
to combine the replica method with infrared perturbation theory to compute the bipartite entanglement entropy of a 
segment of length 2L with the rest of the system in the case of a CFT perturbed at its boundary by the interaction 
with an impurity, namely the Interacting Resonant Level Model, or the one-channel anisotropic s = 1/2 Kondo model. 
Our approach allows for a determination of the Renyi and entanglement entropies in a systematic expansion around 
the strong coupling fixed point that corresponds to the limit L>T B , where T B is the boundary energy scale generated 
by the interaction with the impurity. We have also investigated the case of a finite size system, and shown that a 
similar expansion can be carried out in this case. 

Our results indicate that generically, the effect of the interactions on the entanglement entropy is weak: the first 
contribution to the entanglement entropy that depends on the interactions scale as 1/(LT B ) 6 at infinite size, and as 
1 / (LqLT^) at finite size. On the other hand, we emphasize that the effect of interactions is strongly enhanced in the 
case of strong coulombic repulsion U > in the one-wire IRLM. 

Note that our method could also be extended to more complicated geometries, like finite temperature and finite 
size, the only (but significant regarding the difficulty in computing the correlators!) difference being that the orbifold 
theory should be put on a torus. 

At the end of such a technical paper, it is of course natural to ask how the perturbative expansion compares with 
say a numerical evaluation of the quantity S lmp over the whole flow from UV to IR, and how, in particular, the UV 
limit - which is non perturbative in 7 - is approached. This will be discussed elsewhere [18]. Similarly, it is clear 
that in most circumstances, the quantity of physical interest will not be be S lrap itself, but rather the entanglement 
between the two wires connected by the resonant level. This will involve more technicalities, and will also be discussed 
elsewhere. 

Acknowledgments: We thank Ian Affleck, John Cardy, Pasquale Calabrese and Peter Schmitteckert for discussions. 
HS was supported by the French Agcnce Nationale pour la Recherche (ANR Projet 2010 Blanc SIMI 4 : DIME) and 
the US Department of Energy (grant number DE-FG03-01ER45908). LF acknowledges the support of the French 
Direction Generate pour l'Armement. 



Appendix A: Twist fields and orbifold representation 



We consider a CFT of central charge c (later one will apply this to the case c = 1) with fields living on the surface 
TZn with branching conditions between the ./V sheets of the surface. As argued in Refs. [31, 33], this is equivalent to 
considering a theory leaving on the complex plane but with twisted vacuums |orb) = |0) and (orb| = (0| $(u) 
and central charge c = Nc. The equivalence as the following meaning: given a local operator X(z^') in the original 
theory defined on TZn (here z belongs to the punctured complex plane C\{u, v} and j = 1...N is a sheet index), 
one can form linear combinations XW(z) = £\ e 2mjk / N X(z^) that carry a well defined charge under the Zjv cyclic 

exchange of the sheets. The equivalence states that to each field X^ there corresponds a field X^ in the twisted 
theory such that the correlators of X^ fields on TZn coincide with those of X^ according to: 

(...*(%)...) = <^&h^MM (A1) 

\ 'n N ($(«)$(«)) 

and in the following we will denote X^=X^ this equivalence between operators. 

The twist (and anti- twist) fields $ (and $), are primary fields of conformal dimension h = -!k ^ ■ (Note that in 
this paper the CFT one considers are chiral ones as the result of our unfolding of the original problem, so that h is 
not defined). |0) is the vacuum of the chiral theory living on a single sheet. The two point function on the plane of 
these twist fields reads : 

'$(«)$(«)) = l/(«- v) 2h (A2) 



Appendix B: Energy momentum tensor in the orbifold theory 



In this appendix one derives some identities for the operator product expansion of the energy-momentum tensor 
and its " orbifold powers" , by using the equivalence between operators in the untwisted sector of the orbifold theory 
and Zjy-neutral combinations of operators in the original theory living on the Riemann surface TZn- 
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1. The orbifold energy-momentum tensor at infinite size 

We start by considering the Zjv-neutral combination T^{z) = J2f=i Tn N (z^) = Y^j=i Tj(z) 01 the stress energy 
tensor (here j labels the sheets of the Riemann surface Hn), and show, following Ref. 30, that it coincides with the 
orbifold energy- momentum tensor T or b (and we will write TLrb^X^ ^i) m the sense that 



A orb $(«)$(«) 

= V , ~ v ' (131) 
$(u)*(w), 



where X = T(°)(^i)...T(°) (z fe ) is a collection of energy-momentum tensors inserted at arbitrary points on <C\{u,v}. 
The left hand side of Eq. (Bl) (first route) can be computed by using the uniformizing map that sends TZn onto the 
punctured complex plane: 



z^e = >u = e> ■ (B2) 

z — u 

On the other hand (second route), the right-hand side of Eq. (Bl) can be computed using a Ward identity involving 
the energy-momentum tensor T or b- 



a. One point correlator 



Let us follow the first route to compute the one-point function of ouTZn- Under the uniformizing transformation 
(B2), the energy-momentum tensor transforms as: 



T nN {z) = (d z co) 2 T c {u) + — 2 



uj- 



so that 



(%W> = « 



d z <jj 
to 



24 v 1 



a-o 4 



[N(v-u)] 2 i 2 



(B3) 



(B4) 



where one used d z uj 



N ^ v _ u ^ — f — • Since this one-point function doesn't depend on the Riemann subsheet where 
it is evaluated, one thus gets the final result by just multiplying by N: 



T<®(z] 



(u — v) z 



(B5) 



N(u - v) 2 i 2 N [z- u) 2 {z - v) 2 ' 
According to the second route, one also has: (T^(z))„ = w i^ u t , the energy-momentum tensor 

* ' /vjv (&(u)&(v)) ' 



of the orbifold theory with central charge c = cN , and ^$(m)$(u)^ = (u — v) 2h where h is the conformal dimension 
of the twist fields. As noted in Ref. 30, one can indeed recover the result (B5) by using the fundamental OPE: 



M>(u) d u ®(u) 
T olh {z) .$(«) = + -2-U. + reg. 



(z - u) 2 



that leads to the following Ward identity: 
[T olh (zMuMv)) 



z—u z—v 



(z-u) 2 (z— u) 2 

where the conformal dimension of the twist h is identified as: 

h = k/N = — ( N - — 

1 1/1 \ AT 



$(w)$(v) 



h(u — v) 



2- 2h 



(z — u) 2 (z — v) 2 
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N 



(B6) 



(B7) 



(B8) 
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Note that we can now easily compute the simplest contribution — J _ Q dz (T(z)) . First one notices that 



in this expression, T can be replaced by T since the term T — T oc dQn N is a total derivative and immediately 

h (u— v) 2 Ainh 



A {0} 



dz 

q N (z — u) 2 (z — v) 2 u — v 



(B9) 



where one recalls that u = —iL and v = iL, and where a factor of N in the last expression comes from the N 
disconnected part of the path of integration. An easy way to compute the integral in (B9) consists in deforming the 
contour of integration in the half plane Im(z) > and pick up a residue at z = v. 



b. Higher order correlators 



One can explicitly check the equivalence of the two routes when considering higher order correlators. Let us for 
example compute the two-point function of T or b ; one gets 



£< r ii(i) r * = E 



31 ,J2 



31 ,32 



(i-ei) 4 q~6) 4 
(n{u-v)y (ae 2 ) 2 



c [oj^w^f 



l (i-6) 4 (i-6) 4 



N 2 (66) 2 («-«) 4 
where one has written ui^ — e 2m ^' N oj a , and one has introduced the function 



* 2 + 7j/ 2 ,2(^M) 



(BIO) 
(Bll) 



fp,q( X ) " ^ E (1 



e « 



■(32- jl) 



3l,32 



(1 _ xe -ff(32~3l)y+q 



(p+q-l)\ y \y N -l) 



y=x 



(B12) 



The last equality shows that /„„ is in fact an entire function of x N ; it can be derived by using an integral representation 



fpq(x) = jijr j c ± ydft-i) (i-^)p+g wnere circles the unit circle, clockwise at radius 1 — + and counter-clockwise 



N(N 2 -1) 



Note 



at radius 1 + + . Moving the contour one gets Eq. (B12), so that f2,2(x) = N 3 ( ^J^at^ ) H g (i-x N ) 2 ' 

that, as it should, the final result (Bll) does not depend any longer on the sheet the points z\ and z 2 are chosen in. 
Putting everything together, one gets: 



£ (T n { Zl )T 32 {z 2 )) nN =h 2 D 2 



cN 



2hD 



(z 12 ) 2 2( Zl2 y 



with D = 



(v-u) 2 



(ZX—U)(Z1-V)(Z2-U)(Z2—V) 



(B13) 



This is again consistant with the Ward identity 

fr orb (z 1 )T orb (z 2 )$( u )$( u ) 
+ 



2(zu) 

h , h 



4 (*(«)*(») 



2 , 9,i 



r orb (2 2 )$( u )$( w ) 



(B14) 



2. Computations at finite size 

As we have seen in the section V, in order to compute the correlators on Cn, with the equivalence with the orbifold 
theory, one has to make another change of variables from the plane to the cylinder z — ^rlnu;, that change the 
integration contour and thus forbid to move it and wrap it around the twist fields, as we did in the infinite size case. 
For instance 



E 

i=l 



T(z«) ) = X V ; „ ; V " c (B15) 



($(tt)$(u)) 

((g- 2 r orb M$(ii)i(?)) c g /& 
<a>(S)$(^)) c 11 ' 



12 V 



{z-M 
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And without too much pain, one can manage to compute the very first terms of the contributions to the Renyi 
entropy, up to order 0(T~ 2 ) (remember that <9<p( _fe ) = dip^* and v — e %7rL ^ L ° — u*) 



A 



7T 

Lq Jo 



hu;(u - 



\u-v\ 



N 
2M> 



2nh 



fi{L/L ) 



dco- 



24L w 



div. 



(B16) 



where 



A {0;0} = /AO} 



du\duj2 



hir 2 LJiU!2 



(u — vY 



Nir 2 uiiUJ2 



Ll(ui - cu 2 ) 2 |wi - v\ 2 \lo 2 - v\ 2 4Lg(wi - w 2 ) 4 



2h7r 2 f 2 (L/L )/L 2 



div. 



v- 



2h {2hf 1 (L/L ) 2 + f 2 (L/L )} + div. 

/i(a) = a (1 + 7r(l — a) cot wa) 
/ 2 (a) = /i(a) 2 + a 2 (l-a) 2 7r 2 



(B17) 

(B18) 
(B19) 



Note that these integrals are well defined, even in the presence of the (coj — Uj)~ k terms, because the integration 
points are " avoiding each other" . This is a consequence of our point-splitting regularization scheme. 

As discussed in the main text, the diverging terms in the integrals stemming from "intricated" integrands are 
proportional to N only and are cancelled out by terms stemming from the denominator of Eq. (20). 



Appendix C: The composite field 04 



In this appendix we establish the operator product expansions of the first composite field in the neutral sector 
of the orbifold theory with the twist operator $. The orbifold form of the perturbing operator is: 



04 Orb M=E( T ^)M : f j = T j + adQ j 



(CI) 



3=1 



where z is a coordinate on the complex plane (punctured at points where the twist fields is inserted) , and on the right 
hand side the index j indicates the sheet of 1Z n- The parenthesis in the second line of (CI) indicates normal ordering 
(A.B)(z) = 2^ $ ^fA(z + x)B(z) that is performed on TZ N . 

Introducing the Virasoro generators' action on a local field 0{z) via L-k{z)0{z) = =i- <£_ dx (x — z) 1_fc T OI ^(x)0(z) 
[34], one parametrizes the OPE of 0° rb and $ as follows (only the descendents of the twist field can appear): 



0? b (z) • $(0) 



a.iL _ 1 +b 2 L-2 a 1 L_ 1 +b 1 L- 1 L-2+ciL- 3 



*(0) 



(C2) 



We will determine the coefficients ai,bi,Ci (see Eqs.(C14)) following the simplest route: one starts by com- 
puting correlators Gxo = {X{zi)0{z)) nN where X{zi) is an insertion of a collection of stress energy tensors 
T or b at points Zi, and one then expands this expression when z — > u. On the other hand, one has Gxo = 

(^$>(u)$(v)\ (^X(zi)G(z)$(u)$>(v)^: by identification this fixes constraints on the OPE C(z)-$(u) that are sufficient 

to determine the a priori unknown coefficients in Eq. (C2). 
The operator (CI) can be split in 3 pieces, 

0? h =(T 2 ) 0lh + a ({dQ, T}) orb + a 2 {dQdQ) mh 

Let's start with (T 2 )orb(z)=Y,j( T j T j)( z ) = § z J2j T j( z i) T j( z )- Notc that the summation implies that (T 2 ) orb ^ 
(Torb) 2 - Starting from the correlator (BIO), inserting 5{j\ — 32} in the sum, and performing the normal order one 
finds: 



N 3 (u- 



- 2 (i-0 2 



GO 



Pi(H) 



(C3) 



22 



with Pi(X) = (1 - X) 2 (29N 2 - 11) + 9CW 2 (1 + X) 2 

Parametrizing the OPE of (T 2 ) or b with the twist field as: 



(T 2 ) OTh (z)-m 



a 3 L_i . a 2 L_ 1 -\-b 2 L_ 2 . a 1 L'_ 1 +b 1 L_ 1 L„ 2 +c 1 L„ 3 



*(0), 



(C4) 



one can then expand the correlator (^(T 2 ) or b(z)<fr(u)$(v)J when z approaches one of the twist fields, say <&(«), and 
identify the two expansions obtained via the OPE (C4) and via (C3). The first expansion involves correlators: 



$(u) L_ii_ 2 $(v)\ = 



dz /„ , ,~. . \ (n + 
— (T olh (z)®(u)$(vy 



2iir J v (z — v) 
6h(h + 1) 

( u _ w )3+2/> 

Now expanding (C3) when z — > it, one gets the following identifications: 

1 

5o/V 
1 



(u — v) 



n+2h 



(n > 0) 



(C5) 



«4 

as 



[fc(22 + 5c) + 9cN] 
[2/i(22 + 5c) + 3ciV] 



2(2/i + l)a 2 + 36 2 
(2/i + l)(2h + 2)ai + (3/i + l)6i + 2ci 



5ciV 
^(22 + 5c) 

^(22 + 5c) 



(C6) 



In order to remove the ambiguity in the determination of the unknown coefficients aj,6j,Cj, one needs to consider 
higher order correlator with more insertion of T^: 



E (Tj{zi){T k T k )(z 2 )) U[ 



(i-ft) 4 (i-6Q* 

N*(u- V )ni$ 



E 

3,k 



1440(1 - £ 2 ) : 



■p 2 (6 



KC Ulfui 2 



U> 



12 



2 9 4— p 

i2(i-6) 2 p ^ W42; 



(i-a) 4 a-6) 



8 



N\u - vf Q$ 



E 

3,k 



1440(1 - £ 2 ) 



+ KC / 2 ,2 — 



with 



P 2 (X) = (-22 - 100iV 2 + 122^ 4 )X + (ll-130iV 2 + 119^ 4 )(l+X 2 ) 

Q {X) = 36(1 -X) 2 

Qx{X) = -m{l-X)[{N+l)X + N - 1] 

Q 2 pO = 5(l-X) 2 + 187V(x 2 -l)+iV 2 (13X 2 + 10X + 13) 



(C7) 



(C8) 



This fixes the coefficients a 2 ,6 2 . A last constraint is needed: we compute the 3 points correlator 
J2j k i (Tj(zi)Tk(z2)(TiTi)(z3)) nN (its expression is not particularly enlightening so we do not display it) and we 
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take the limit z 3 — ¥ v. Finally we find: 
22 + 5c 



23ch - 22h 2 



a i 



Cl 

«2 

a?, 

«4 



5c (2 + c - 7h + ch + 3h 2 )(c - lOh + 2ch + I6/1 2 ) 
22 + 5c 2(c 2 - 19c/i + 3c 2 h + 54h 2 - 2\ch 2 + 2c 2 h 2 - 90h 3 + 22ch 3 + 48/i 4 ) 
~~bc (2 + c - 7h + ch + 3h 2 )(c - lOh + 2ch + 16h 2 ) 

22 + 5c 2(h - l)(c - 12h + ch + 3h 2 ) 

5£ 2 + c - 7h + ch + 3h 2 

22 + 5c 2h(c + 8h) 



5c c - lOh + 2ch + I6/1 2 
22 + 5c 2h(c - 22h + 2ch + I6/1 2 ) 



5c 

3 „, 22 + 5c 

l + 2k — 
9 , 22 + 5c 

5 +k — 



Wh + 2ch+16h 2 



(C9) 



So far, the computation was totally general and applies to an arbitrary CFT put on a In ■ In our case, one deals with 
a free boson Q cx id(f>+, so that one should put c = 1 and c = N in Eqs.(C9). 

We are not done so far: there is a second piece to 0^ h ', namely terms that depend on the deformation parameter 
a. On the plane, this deformation parameter simply amounts to a change of the central charge: T is an energy- 
momentum tensor of a c = 1 — 6a 2 CFT. However, it would be wrong to identify Of b with (T 2 ) orb , the orbifold 
version of the energy-momentum tensor of the deformed theory. Instead, what one should do is put the c = 1 theory 
on the orbifold, and determine the orbifold version of the operators a 2 J2j (dQjdQj) and q^] - ({dQj, Tj}). 



Consider ■ (dQjdQj): the same game can be played with it, its OPE reading: 



(dQdQ) olh (z) ■ $(0) 



3L-1 



a,2L_ 1 +b 2 L-2 a 1 L'_ 1 +biL-iL-2+ciL-3 



$(0) 



(CIO) 



Proceeding in a similar way (we do not reproduce the steps, that are identical), we obtain the coefficients: (a,, bi, Cj) 



^(ai^ijCi), except for a 3 = -ga 3 



we find the remarquable relation: 



I and 04 = — |a4 + 2. Comparing this with the OPE: 



E 



(dQjdQj 



(Wj) 



(Cll) 



(C12) 



so that keeping only terms that are even in a (terms linear in a are easily showed to be total derivatives, ({T, dQ})Tz r 

2 



^d(TQ)n N , and in general are odd under Q — » — Q, and they don't contribute to expectation values): 



(9° 



a — even 



= ( 1 - ) ( T ")orb + ^d 2 T mh 



(C13) 



Restricting again to terms that are even in a, one gets: 

{Si,bi,Ci}i<2 = {\a,i,\bi,\ci} 
a 3 = Xa 3 + 2a 2 /3 
S4 = A 04 + 2a 2 

with A = 1 — 4a 2 /9, and {dj, 6j, c^} are given by Eq.(C9) with c = 1. 



(C14) 



Appendix D: Computations with a direct quantization 



We hereunder report some details useful in order to understand how to explicitly compute correlators with the 
mode decomposition. 
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1. Energy density on the plane 

With the mode decomposion it is easy to show that (tp^ (£)d^~ fc ) (£')) = Y^ m - c "F 7 ^ I — 

g (k) (LO = (dv {k) (tW- k) (0) 

_k_ 

=- r ( i-* )+ g* 

4tt(£-£') 2 L N' £ N 
we then have the expression for the /c t/l — energy-momentum tensor (T 1 ^^)) — —2n (g^(^ + e,^) + 



so that we find 



4iV 2 £ 2 

and summing up this expression over k to get the total energy-momentum tensor, we have 

c(N 2 - 1) 1 



(T(0) = 



24iV 2 £ 2 



(Dl) 
(D2) 

(D3) 
(D4) 



with c = iVc (c — 1), since initially we have N real bosonic fields. This is the result that we would get doing the 
conformal transformation uj = , this expression is the value of the schwartzian derivative. The branch point in 
£ = induces stress that leads to a non vanishing expectation value for the energy-momentum tensor. 



2. Cut of finite size on the plane and on the cylinder 



If we are interested in a cut between two points u, v on the z— plane, all we have to do is make a conformal 
transformation that sends v — > and u — ¥ oo. Namely 

{ - — (D5) 

z — u 

since dtp^ is a primary field, it is easy to compute with the correlator (D2) 

and summing up over k, we find the expectation value of Eq. (B17) 

{{)) 24iV2 {z-uf{z-vY [ ' 

we see in the next section that we can find directly this result using the energy momentum tensor's transformation 
identity for since the schwartzian of the transformation (D5) is equal to zero. 

Now we want to compute correlations on Cn ■ cylinders of size 2Lq linked by a finite cut between u — —iL and 
v — iL. To come back to the usual situation where the fields acquire a phase when turning around z = 0, we have 
to make the following transformations : first explode the cylinder on the plane, where u = eT mL l L ° , v = e l7rL / L ° and 
then transform v — !• and u —¥ oo, namely it is given by 

p TTz/Lo _ +iirL/L 
Z = —^7 ^TTT (D8) 



where z £ C. 
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3. Schwartzian 



It is interesting to check explicitly that even in this twisted geometry, at hnite size, and with this mode decomposition 
one finds that the schwartzian appears when computing 



T^(z)) =lim-2vr 

/ cyl. e^O 



d^ k) {z + e)d^- k) {z) 



= -2tt 



dz 



{u,v} I 



{0,oo} 1 



(z)f z (z + e)(dv (k) (az + t)W~ k \az))) / + 



n 2 k (1 - k) e 27rz / L ° (u-v) 2 or 2 
4:N 2 L 2 (e* z / L ° - u) 2 {e* z / L « - v) 2 ~ 24L 2 



(D9) 



so that we will now directly use the well-known identity : T' (z) 
momentum tensors. 



T (0 + ilU; z}i even for the energy- 



4. Example with (T 2 ) 

Where the quantization scheme becomes really useful is when one is interested in computing expectation values of 
operators that have no simple expression in their orbifold version, for instance (T 2 ). 

All we have to do is to use the mode decomposition 56 of the basic fields tpu\ and use the Wick theorem to compute 
explicitly 

f dz((f 2 )(z)) =W dz<(2?)(*))£; u} +a 2 X; [dziidQjdQAW^ (D10) 
J x =o x ' c « ; j x =o ' ~r J 



where only the non vanishing expectation values have been written. The integrand can be expressed with the fields 

d(p^ as 



+ 2(-2tt) 2 (d^ k \z + e) (d<pM(z)d(p(- k ">(z + e)) 



(Dll) 



Term A The term A expressed in the w— plane (change of variables 46) is given by (we denote a = L/Lq) is 

h 2 ir 3 (u- v) 4 uj 3 hit 3 {v-u) 2 LO Nn 3 



A = du 



L 3 N 3 \d-u\ & 12L 3 N\u - w| 4 576L 3 uj 

div. 

h 2 n 3 a 2 / „ A (a) „ , , , . \ /i7r 3 2 „ . . 

' 15^4H- - 6/i(a) - 5a ) + — a 2 h{a) + div. 



3NL 3 V sin 2 



(D12) 



Term B The term B gives 



B 



B = 



\22h 2 ir 3 (v-u) i uj 3 hTT 3 {u-v) 2 uj ( 25 , 90w(w + w) ( 54(2 -vf 



5NL 3 \uj 



30L 3 



lliVTT 3 

1U0L 3 lj 



div. 



22h 2 n 3 a 2 ( /i(a) 



15AL 3 V sin 2 



15- 



6/i(a) - 5a 



15L 3 



a 2 /i(a) + diu. 



(D13) 
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Term C And finally, the term C gives 

12h 2 ir 3 (u- v) A u 3 hn 3 (u-v) 2 uj 



C 



Ah 2 n 3 a 2 
5NL 3 



5NL 3 \uj 



5L 3 



10uj(u + v) 6(u~v) 2 lj 2 



Nir 3 
240L 3 oj 



div. 



15 



/i(«) 
sin 2 air 



6/1 (a) - 5a 



2hn 3 
5L 3 



a 2 /i(a) + div. 



(D14) 



The IR divergences that come when integrating at large distance from the insertion of the twist operators will 
actually not depend on the twisted geometry. Hence these divergences will only get multiplied by a factor N when 
one integrates on TZn, so that they will be compensated at every orders. This is verified for the diverging terms 
stemming from A and B which are indeed linear in N, so that the expansion in power series of 0(1/LT B ) of the 

exponential exp(— efeA^J) will be exactly compensated from terms of the denominator exp(— N dzA^) in 
the expression of Rn- 

Combining the results one finally gets (a = Lq/ L) 



a{1} = (9 - 4a 2 ) 



3h 2 n hK 3 a 2 h{a) 

— Mia) H 

NL 3 ' v ' 



2L 3 



(D15) 



where 



/ 4 («) 



15 



(l5/i(o)/sin 2 (a7r)- 6/1(0) -5a) 



= 1 - 27r 2 a 2 /5 + 7r 4 a 4 /15 + d(a 6 ) 
It is easy to check the result known at infinite size for AW (from Eq. (40)): 



{1} 3(9-4a 2 )fe 2 vr 
5A^ 3 



3vr/i 2 (4-L>)(l-4L>) 
DNL 3 



(D16) 
(D17) 

(D18) 



Note that now, gives non- vanishing corrections to the entropy (terms proportional to h). These corrections 
are due to finite size effects, they are proportional to a 2 and depend on the interactions (a 2 term). 
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